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ABSTRACT 

We present kinematic measurements of the thick and thin disks in two edge-on galaxies. We have 
derived stellar rotation curves at and above the galaxies' midplanes using Ca n triplet features mea- 
sured with the GMOS spectrograph on Gemini North. In one galaxy, FGC 1415, the kinematics above 
the plane show clear rotation that lags that of the midplane by ^20-50%, similar to the behavior seen 
in the Milky Way. However, the kinematics of the second galaxy, FGC 227, are quite different. The 
rotation above the plane is extremely slow, showing <25% of the rotation speed of the stars at the 
midplane. We decompose the observed rotation curves into a superposition of thick and thin disk 
kinematics, using 2-dimensional fits to the galaxy images to determine the fraction of thick disk stars 
at each position. We find that the thick disk of FGC 1415 rotates at 30-40% of the rotation speed 
of the thin disk. In contrast, the thick disk of FGC 227 is very likely counter-rotating, if it is ro- 
tating at all. These observations are consistent with the velocity dispersion profiles we measure for 
each galaxy. The detection of counter-rotating thick disks conclusively rules out models where the 
thick disk forms either during monolithic collapse or from vertical heating of a previous thin disk. 
Instead, the data strongly support models where the thick disk forms from direct accretion of stars 
from infalling satellites. 

Subject headings: galaxies: kinematics and dynamics — galaxies: formation — galaxies: structure — 
galaxies: spiral 



1. INTRODUCTION 

Sinc e the discovery of the Milky Way's thick disk in 
1983 ( Gilmo^^OIiIdl[l983), there has been a growing 
body of evidence that the thick disk is a structurally, 
chemically, and kinematically distinct component of the 
Galaxy. Structurally, the Milky Way's thick disk has a 
significantly larger scale height than th e thin disk (for re- 
views of the thick disk scale h eight se e lReid fc MaiewskH 
(|l993); iBuser et alJ ((1999^1 : iNorris! l)1999|) and refer- 
ences t herein), and a possibly somewhat longer scale 
length llRobin et aUllQM lOihall200lT: iChen et al1l200ll 
lLarsen fc Humphreys! I2003|) . Chemically, thick disk 
stars are more metal-poor and older than stars in th e 
thin disk lIReid fc Maiewskill7T)Q3l ICIfiha fc Beersll2000ft . 
They are also significantly enhanced in a-elements, 
compar ed to thin disk stars o f comparable iron abun- 
dance l |Prochaska et^i^ [2QQ|j: ^LUt^isign^^L| 2001; 
i Bensbv et al.ll2003UFeltzing et alJl2003HMishenina et all 
l2004HBrewer fc Carnevll2004|) . Kinematically. thick disk 



stars have both a larger velocity dispersi on and slower 
net rotation than stars in the thin disk (Nisscn li 



Chiba & Beers 2000; Gilmo re et al.l2 002: Sou biran et al.l 
120031 iParker et alJ!2004|) . All of these facts lead to the 
conclusion that the thick disk is a relic of the young 
Galactic disk. As such, it provides an excellent probe 
of models of disk galaxy f ormation (see recent reviews 
bv iNissen et alJ pOO^ : [E eeman fc Bland-Hawthornl 
(2003)). 

Originally de tected in SO galaxies l)Bursteml 119791: 
iTsikoudil Il979|) . thick disk components have since 
been detected in many galaxies, including SO's 
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ide Crjis fc van der Kruil H£9J Ide Griis fc PeletieH 
19971 IPohlen et al.l I2004h . Sh's <lva,n der Krnitl ' 



Shaw fc Gilmore! 119891 Ivan Dokkum et al 
iMorrison et all Jl997t IWu et all J2002D. and later 
type galaxies l| Abe et a l.l Il999t iNeeser et all 120021: 
lYoachim fc Dalcanton! 12005^ . Across all Hubble types, 
photometric decompositions consistently find that thick 
disk scale heights are 2-6 times larger than thin disk 
scale heights and t hat the thick di s ks' sc a le lengths 
are sl ig htly larger |Shaw fc Gil more! (11 989th IWu et alJ 



Neeser et all (120021) : IPohlen et alJ 
Abe et a.l.ll|1999fl ) 



(2004), 



i200: 

see also 

Thick and thin disks are likely to have distinct for- 
mation mechanisms, given the systematic differences be- 
tween their properties. The structure, dynamics, and 
chemical abundance of the thin disk strongly suggest that 
the majority of its stars formed gradua lly from a rotating, 
disk of high angular momentum gas llFal l fc Efst athiovJ 
ll980HMatteucci fc Francoisll989tlChiaDnini et alll997fl . 
In contrast, the formation of the thick disk is still poorly 
understood and is likely to be more complex. The large 
scale height of the thick disk suggests that its stars 
were either (1 ) vertically "heated " from a previously 



1999: iRobin et all 119961 IChen et all 
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from gas with a large scale height I 




1962; 


Norris & Ryan 1991; Burkcrt ct al. 1992: Krouoa 


2002; 


Fuhrmannl |200# IGilmore & Wvse! I1986D. or (3) di- 



rectly deposited at large scale hei ghts during the accre- 
tion of smalle r satel li te galaxies llBekki fc Chibal 2001; 
Gilmore et al.l 120021: lAbadi et all 120031: [Martin et alJ 
2004HNavarro et al.l2004HBrook et al.l2004|) . Most mod- 



els of thick disk formation fall into one of these three 
cases. In the first, the thick disk stars form initially in 
a thin disk. In the second, they form within the thick 
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disk itself. In the third, they form outside of the Galaxy 
entirely. 

All three of the above scenarios are compatible with 
the structure of the Milky Way's thick disk. However, 
they all imply very different behavior for the early evo- 
lution of disk galaxies. The first case would suggest that 
disks form largely through smooth gas accretion that is 
occasionally punctuated by minor merging events. The 
second case would suggest either that disks form pri- 
marily through smooth monolithic collapse, with thick 
disk stars precipit ating out of the collapsing gas cloud 
l)Eggen et al]ll96 2'). or that the thick disk forms from a 
disk of ga s that has bee n energetically heated by star 
formation l)Kroupall2002j) . The third case would suggest 
that disks can form entirely from merging sub-units, in 
spite of their highly ordered present day structure. The 
degeneracies between these models must be broken be- 
fore the properties of the thick disk can be used as a 
constraint on theories of galaxy formation. 

Past attempts to distinguish among formation models 
have often relied on the chemical abundances of thick 
disk stars. The low metallicities of thick disk stars 
suggest that they formed early in the evolution of the 
Galaxy. Their a-enhancement suggests that the time 
scale for their formation was sufficiently rapid to suppress 
iron enrichment by Type la supernovae. Unfortunately, 
the resulting implication that the thick disk formed early 
and rapidly can be accommodated by all of the possible 
formation scenarios. iBland-Hawthorn fc Freeman! ( 2004) 
suggest that metallicity measurements of more than 10 5 
stars would be required to definitively constrain the ori- 
gin of the Milky Way's thick disk. Metallicities of thick 
disk stars therefore remain a somewhat weak discrimi- 
nant among competing models. 

In contrast to metallicities, stellar kinematics show 
more promise for constraining the origin of the thick 
disk. If the thick disk results from dynamical heating 
of a young previously thin disk, both components should 
have similar rotation curves, assuming similar angular 
momenta for the gas that falls in before and after the 
heating event. If the thick disk forms from gas as it col- 
lapses into the galaxy, then angular momentum conser- 
vation requires that the rotation of the thick disk must 
always lag the rotation of the thin disk, by an amount 
that depends on the relative scale lengths of the thick 
and thin disks. Finally, if the thick disk is assembled 
from stars originally formed in merging satellites, then 
the relative kinematics of the thick and thin disks are 
likely to be highly variable. 

The rotation of the Milky Way thi ck disk is marginally 
consistent with all of these scenarios jBjndce^e^alJl^S 
Robin et al Jl993 iBuser et alJll999UGilmore et al.ll2002t 
SoTabTnm'etafTl^OOffi . Measurements have consistently 
shown the Galactic thick disk to be co-rotating with the 
thin disk, but with modest lag. The exact value of the lag 
varies by a factor of two depending on the study, rang ing 
from - 40 km s " 1 (iR.eidl 119981 iChiba fc Beersll2000ft up 
to ~ 90 km s" 1 l|Gihnore et al.ll2002HParker et al.ll 20041. 
The smaller rotational lags can be easily produced in 
most thick disk formation models, but the larger lags 
observed recently are more difficult to explain without 
satellite accretion. However, rather than ruling out out 
alternative models, the larger lags are often interpreted 
as contamination from a single disrupted satellite. 



Although data within the Galaxy have not been suf- 
ficient to definitively constrain the origin of the thick 
disk, the Milky Way is only one single case. Instead, 
measurements of thick disk kinematics in several galaxies 
are necessary to exploit the full power of kinematic dis- 
criminants among formation models. Scenarios involving 
mergers will naturally lead to much more kinematic vari- 
ation among the population of thick disks than if thick 
disks formed during monolithic collapse. Even among 
merging models, we would expect the relative kinemat- 
ics of the thick and thin disks to vary more if the thick 
disk stars are directly accreted than if the thick disk was 
built from vertically heated thin disk stars. By measur- 
ing the relative kinematics of thick and thin disks in a 
large sample of galaxies, one can therefore discriminate 
among thick disk formation models far more effectively 
than by studying the Milky Way alone. 

In this paper we present the first results of our 
on-going study of the kinematic properties of thick 
disks in a large sample of late-ty pe edge-on galaxies 
ijDalcanton fc Bern stein 2000, 2002). We have measured 
of the kinematics of thick and thin disk stars in two edge- 
on late-type disk galaxies. We first present the sample 
selection and observations in We then describe the 
data reduction and extraction of the rotation curves in 
521 In we decouple the observed rotation curves to 
derive the rotation of the thick and thin disks. In 
we measure radially resolved velocity dispersion profiles 
above and below the plane. We then discuss the im- 
plication of our initial results for models of disk galaxy 
formation in i^and summarize our conclusions in §7\ 

2. OBSERVATIONS 

We have carried out long-slit spectroscopic 
observations of two galaxi es drawn from the 
iDalcanton fc Bernstein! (J2000) sample of edge- 
on late-type galaxies. The original sample of 49 
galaxies was selected fro m the Flat Galaxy Catalog 
( Karachejrjsej_el_aDll22 3l) and i maged in B, R, and K x 
I Dalc anton fc Bernstein! I2000fl . IDalcanton fc Bernsteinl 
(2002) used this imaging to demonstrate the ubiquity 
of thick disks around late-type galaxies. We have since 
used two-dimensional decompositions of the galaxy 
images to measure structural parameters for the thick 
and thin disks. This deco mposition will be presented i n 
detail in a future paper l)Yoachim fc Dalcanton! i2005) . 
and we now give a brief overview of our decompositions. 

When fit ting models t o our galaxies, we adopt the 
method of iKregel et~aT] (|2002) and use a Levenberg- 
Marquardt % 2 minimization routine. Each disk is mod- 
eled as having a surface brightness 

z) = ^ {R/h R )K 1 {R/h R )scch 2 (z/(2h z )) (1) 

where K\ is a modified Bessel function of the first order, 
So is the edge-on central surface brightness, R is the pro- 
jected radius along the major axis, and h z is the expo- 
nential scale height. We weight individual pixels by the 
inverse model. This method allows one to fit regions of 
very low surface brightness without being overly skewed 
by brighter midplane features. Comparing the formal re- 
duced x 2 values of the single dis k and two-disk models 
using an i^-test (Matthews 2000), we find that the two- 
disk models are strongly justified. To derive estimates 
for our parameter uncertainties, we fit a series of models 
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Fig. 1. — Vertical residual profiles for photometry fits to our 
galaxies. Plots were made by collapsing residual images along the 
major axis and normalizing by the model galaxy and averaging 
both sides of the galaxy. In both galaxies, the single disk pro- 
vides a poor fit and leaves excess flux at high galactic latitudes. 
Changing the functional form of the single disk does not alleviate 
this problem. We have also plotted the residuals from our extreme 
models. As expected, the "faint thick" model leaves excess flux at 
large z, while the "bright thick" over-subtracts at high z. 



varying the functional form of the vertical profiles (e.g. 
sech(z), sech 2 (z)) for both the thick and thin disk as well 
as comparing models which were fitted with the midplane 
region masked. We estimate our uncertainties by using 
the full range of convergent values in all the models. Our 
final fitted parameters and uncertainties are listed in Ta- 
ble n] Plots of our residuals collapsed along the major 
axis are shown in Figure ^ Our residuals clearly show 
that subtracting off a single disk leaves excess flux at 
high z while our preferred models fit well across the full 
range of lat itudes. This procedure is disc ussed more ex- 
tensively in lYoachim fc Dalcantonl l)2005(l . 

To define a sample for measuring the thick disk kine- 
matics, we selected several initial candidates for spectro- 
scopic follow-up. We chose to focus our initial efforts 
on galaxies with rotation speeds near V c ~ 100 kins -1 . 
These galaxies are rotating fast enough that their kine- 
matics can be measured with moderate resolution spec- 
troscopy. However, they are also low enough mass that 
they fall below the V c ~ 120 kms" 1 transition w here dust 
lanes become common ( Da lcanton et al.ll2004fl . allowing 
us to probe their thin disk kinematics and structure more 
reliably. From the galaxies with V c ~ 100 km s -1 , we fur- 
ther restricted our sample to the most nearby, spatially 
well-resolved galaxies. The resulting sample was then 
submitted to the NOAO Gemini North observing queue. 



The choice of the specific galaxies targeted was then left 
to the Gemini observers, to maximize the chances that 
one of our targets would be visible at any night during the 
observing semester. In this paper we present results for 
the two galaxies that were selected by the Gemini queue 
- FGC 227 and FGC 1415. The properties of both galax- 
ies are listed in Table 1, details of our observations are 
listed in Tabled and images of the galaxies are shown 
in Figure [3 

We obtained long-slit spectroscopy of the Ca II triplet 
(8498A, 8542A, and 8662A) to measure stellar kinemat- 
ics. We chose two slit locations for each galaxy, oriented 
with the slit parallel the major axis of the galaxy. The 
first slit position was located along the midplane, where 
the light is dominated by thin disk stars. The second slit 
position was located 3-4.5 thin disk scale heights above 
the midplane, where 50 — 80% of the light came from the 
thick di sk, based on the thick a nd th in disk decomposi- 
tions of lYoachim fc Dalcantonl l|2005fl . The locations of 
the slits are superimposed on the galaxy images shown in 
Figure [3 Typical surface brightnesses in the midplane 
and high-latitude slit positions were 21.5 mag arcsec -2 
and 23 mag arcsec -2 , respectively, in the i?-band. 

The galaxies were observed with the Gemini North 
GMOS spectrograph using an 0.5" wide slit and the 
R400-G5305 grating set with a central wavelength of 
~8440A. The CCD was binned by 2 in the spatial direc- 
tion during readout. The resulting spectra have spectral 
resolution of 3.4A (FWHM) near the Ca II triplet region 
as measured by CuAr arc spectra. This resolution allows 
lines to be centroided to ^10 kms -1 . 

Accurate placement of our slit is of great importance 
because we are combining images taken over several 
nights. Our offsets from the midplane are relatively 
small, and an error in placement of order an arcsecond 
could compromise our results. In our remote observing 
program, we instructed the Gemini observing specialists 
to align the slit along the galaxy midplane, after which, 
a pre-programmed offset iterator would shift the slit to 
the appropriate off-plane locations. To place the slit, 
the observing specialist first took an exposure with no 
dispersion element to locate the galaxy midplane. The 
longslit was then moved into the light path and another 
exposure was taken. If the slit did not fall on the galaxy 
midplane, the telescope was offset to the correct position. 
Inspecting the returned acquisition images, we find that 
the observing specialists were able to accurately place 
the slit along the midplane within ~ 0.2 arcseconds (1-2 
GMOS pixels). Because our spectra are dithered along 
the slit with offsets of around 5-10 arcseconds, we have 
checked the overall accuracy of Gemini offset commands 
and found that they are correct to ±1 pixel (0.15 arcsec- 
onds). Overall, we find that the slit was always placed 
within 0.21" of our target position which ensures our tar- 
get regions were always in the 0.5" slit. 

Observations of FGC 227 and FGC 1415 were executed 
during the spring and fall observing semesters of 2003, 
respectively. The image quality was in the 70th-85th per- 
centile, corresponding to typical seeing of ~ 1". Because 
our galaxies are large (~ 1'), as are the offsets between 
the two slit observing positions (> 3 — 5.5"), poor seeing 
has little effect on the resulting spectra. The midplane 
slit positions were observed with three 20 minute expo- 
sures for both galaxies. The thick disk of FGC 1415 
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TABLE 1 

Properties of Targeted Galaxies 
Columns: (1) Galaxy name from the Flat Galaxy Catalog; (2) Measured Heliocentric recessional velocity; (3) Adopted 

DISTANCE FROM KARACHENTSEV ET AL. 12000); (4) 21-CM LINE WIDTH AT 50 PERCENT OF THE PEAK; (5) THIN D ISK EDGE-ON fl-BAND 

centr al surface brightness; (6) Thin disk exponential scale length (Structural parameters taken from Yoachim & Dalcanton 
12005)); (7) Thin disk exponential scale height; (8) Thick disk edge-on central surface brightness; (9) Thick disk scale length; 
(10) Thick disk scale height; (11) Ratio of the thick disk luminosity to the total galaxy luminosity. 



Galaxy 
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FGC 227 
FGC 1415 


5493 
1522 


6261 
2683 


107 
87 


21.5±0.5 
20.9±0.3 


10.0±0.8 
14.8±0.9 


1.0±0.2 
1.2±0.3 


22.8±0.8 
22.2±0.3 


9.3±1.3 
20.0±1.2 


2.0±0.4 
3.2±0.3 


0.50 
0.44 



TABLE 2 
Observing Details 



Galaxy 


RA 


Dec 


Slit PA 


Offset from midplane 


% of Light from 


Integration 




(J2000) 


(J2000) 


(deg E of N) 


(") 


(kpc) 


Thick Disk 


(minutes) 


FGC 227 


02:00:56.7 


19:42:26 


9 






20% 


3x20 


FGC 227 


02:00:56.5 


19:42:25 


9 


3.0 


1.3 


50% 


28x20 


FGC 1415 


12:20:27.7 


01:28:09 


146 






30% 


3x20 


FGC 1415 


12:20:28.0 


01:28:12 


146 


5.4 


1.0 


70% 


46x20 



was observed in 46x20 minute exposures over 9 sepa- 
rate nights, for a total of 15.3 hours. The thick disk of 
FGC 227 was observed in 28x20 minute exposures over 
8 separate nights, for a total of 9.3 hours. The spectra 
were dithered along the slit in a random pattern between 
exposures, with typical shifts of 5". 

3. DATA REDUCTION 

The GMOS CCD's were bias corrected using both bias 
images and fits to 31 columns of overscan. Flat fielding 
was performed with the Facility Calibration Unit (FCU) 
and the twilight sky. FCU flats were interleaved with 
the science exposures approximately every hour to ensure 
complete removal of any variation in the CCD fringing 
pattern. Twilight flats were combined and used to cor- 
rect for uneven illumination of the CCD. 

The spectra were wavelength calibrated using night 
sky emission line s . We combined the sky line data of 
lOsterbrock et all l)1996|) with the IRAF sky line list to 
construct an atlas of lines (or line blends) that would be 
bright and centroidable at our instrumental resolution. 
This procedure yielded ~ 40 lines for chip 2 (containing 
the Ca II triplet) and ~ 30 lines for chip 3. We fit the 
resulting wavelength solution with a 6th degree Cheby- 
chev polynomial. RMS errors in wavelength calibration 
were of order 0.3 A, or ~ 10 km s _1 at the wavelengths of 
interest. 

Analysis of the Ca n triplet absorption line features 
requires accurate sky subtraction due to the high den- 
sity of sky lines at long wavelengths. Each column of the 
wavelength-rectified spectrum was fit with a low order 
(n = 1 — 4) polynomial in regions uncontaminated by the 
galaxy. The resulting polynomial was then subtracted 
from the entire column. Unfortunately, in our deepest 
off-plane integrations, clear errors in sky-subtraction re- 
main. These errors are quite small (typically only 5% 
of the sky), and are not noticible in low signal-to- noise 



images. However, they are a substantial source of uncer- 
tainty in the off-plane spectra, where the surface bright- 
ness of the galaxy is less than 10% of the sky and only 1% 
of the brightest skylines. Measurements of an image of 
the slit suggest that these errors result from small-scale 
(10%) variations in slit width, that lead to systematic 
errors when removing extremely bright sky lines. While 
nod-and-shufnc techniques can greatly reduce this prob- 
lem, our targets are too large to effectively use this ob- 
serving mode. 

After sky subtraction, frames were corrected for at- 
mospheric extinction. No absolute flux calibration was 
applied. The spectra were then Doppler corrected to the 
local standard of rest, spatially aligned, and combined us- 
ing a sigma-clipping algorithm to eliminate cosmic rays. 

Before extracting one-dimensional spectra, we re- 
binned our galaxy spectra to a scale linear in x = In A. 
Such re-binning allows Doppler velocity shifts to be per- 
formed as linear shifts in x. The logarithmic bin size was 
set to 25 km/s in order to preserve the total number of 
pixels. 

Spectra were extracted by summing 30-300 pixels per- 
pendicular to the dispersion direction, corresponding to 
radial bins of 4-40". We varied the extraction bin size 
to ensure that each spectra reaches an average signal-to- 
noise ratio of 10-15 per pixel. In the case of FGC 227's 
shallower off-plane observations, we were forced to bin to 
a lower S/N (~ 5) to obtain reasonable sampling across 
the length of the galaxy. Once extracted, the galaxy 
spectra were normalized by dividing by a low-order poly- 
nomial. Error spectra were computed by analyzing re- 
gions of blank sky on the reduced images, and are clearly 
dominated by the sky-subtraction problems discussed 
above (see Figure |3J). In regions where the sky spectra is 
smooth, our noise spectra closely matches expectations 
from Poisson statistics and read noise. Columns that 
showed excessive systematic residuals caused by skylines 
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Fig. 2. — i?-band images of our galaxies. Solid lines represent the location of the slit jaws for our midplane and off-plane observations. 



(> 5<x above what would be expected from Poisson noise) 
were flagged and masked from later analysis. 

4. TEMPLATE FITTING 

The spectra of galaxies contain detailed informa- 
tion on their stellar kinematics. Methods to extract 
this kinematic information include cross- correlating the 
galaxy spectra with a templa te star ijSimkinl 119741: 
iTonrv fc Davis] Il979t iBottemal ll98Sl). or direct l y fit- 
ting a redsh ifted stellar template |Rix &; White! Il992i 
IKelson et al"]l200rt iBarth et al.ll2002lK Direct fitting has 
the advantage of permitting variable weighting across 
pixels. We therefore adopt the direct fitting method to 
allow us to mask bright skylines that overlap some of 
the Ca II triplet features. We fit the extracted galaxy 
spectra with a redshifted and broadened K-giant stellar 
template in the re stframe region of 848 0-87 00 A, similar 
to the methods of IKelson et al.1 1)20001) and IBarth et al.1 
( 2002). In particular, we construct a model galaxy spec- 
tra M(x) as 

M{x) = {S(x + z)<Z>G(x)}P(x) +C(x) (2) 

where S is the normalized template stellar spectra binned 
logarithmically in wavelength, z is a redshift, G is a Gaus- 
sian broadening function, P is a low-order polynomial, 
C is a flat continuum, and ® denotes convolution. The 
polynomial and flat continuum are used here to correct 
for any errors made when normalizing the galaxy and 
template spectra. The order of the polynomial is kept 
small to ensure that absorption features are not affected. 
While the broadening function can be used to measure 
the line-of-sight stellar velocity dispersion, our large ra- 
dial bins introduce additional broadening into our spec- 
tra. We therefore defer an analysis of the stellar velocity 
dispersion to § and choose to hold the width of the 
Gaussian broadening function constant during the fits. 
We find the best-fit redshift and polynomial coefficients 
simultaneously using a Levenberg-Marquardt routine to 
minimize % 2 : 

X 2 = £(Oi-A4i) x Wi (3) 



where Oi is the extracted galaxy spectra and Wi is the 
weighting derived from the noise spectra with masked 
columns set to zero weight. Uncertainties were calcu- 
lated from the covariance matrix and rescaled so that the 
reduced chi-squared equals unity (i.e., we assume our fits 
are of good quality). This step is mostly conservative in 
that it greatly increases error bars on several points and 
leaves most virtually unchanged. Our rescaling also pre- 
vents any residual unmasked systematic errors from the 
sky subtraction from overly influencing our error spectra 
and fits. An example spectra and template fit can be 
seen in Figure 01 The raw extracted rotation curves are 
presented in Figures 0] and [SJ 

4.1. Fitting the Rotation Curves 

Photometry suggests that our sample galaxies have 
both a thin and a thick disk, both of which contribute sig- 
nificantly to the measured dynamics at each slit position 
(TableEJ • We now attempt to decouple the kinematics of 
these two components. We first introduce a convenient 
analytical equation for the rotation curves: 

V(r) = v + v c - 1 , (4) 

where vq is the recessional velocity at the galaxy's cen- 
ter ro, v c is the asymptotic velocity at the flat part of 
the rotation curve, is the transition radius between 
rising and flat sections of the rotation curve, 7 governs 
the degree of sharpness in the transition zone, and a; is a 
scaled radial param eter equal to r t / (r — rrj). This equa- 
tion is equivalent to ICourteaul l)1997fl 's Model 2 with the 
parameter (3 set to zero. This equation is purely phe- 
nomenological and is flexible enough to fit a wide variety 
of rotation curves. 

We create model rotation curves at each slit position 
as follows. We first adopt independent rotation curves 
for the thick and thin disks, VtMck{r) and Vthin(r), us- 
ing the parameterized rotation curve of equation We 
then assume our observations are some linear combina- 
tion of these two rotation curves, i.e., that each velocity 
measurement is given by 

V(r) = fthickV t hick{r) + (1 - fthick)Vthin(r) (5) 
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Fig. 3. — (Top) Raw extracted galaxy spectra (solid) along with the Poisson noise from the sky spectrum (dotted). (Bottom) Example 
extracted spectra of both midplanes, skyline residuals have been masked for FGC 227 resulting on one of the Ca triplet lines being lost. 
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Fig. 4. — Extracted rotation curves for FGC 227, showing the rotation along the midplane (left) and above the midplane (right). Vertical 
error bars represent the formal uncertainties of the template fitting while horizontal error bars represent the range of radii used in binning. 
Points are placed at the flux-weighted mean bin position. The solid curve is a fit to the data points in the midplane rotation curve, shown 
for reference. Note the very slow (or complete lack of) rotation seen above the plane. 



where fthick = L thlck / L tota i is the fraction of luminosity 
from the thick disk at the slit position. 

We fix the relative contributions of the thick and thin 
disks at each slit position using the photome tric de- 
compositions from lYoachim fc Dalcantonl ((2005T) : as de- 
scribed in S0. We adopt four different possible models, 
shown in Tables [21 & E] and Figure ^ The first is our 
"preferred" model, corresponding to our best fit photo- 
metric decomposition. We also include two extreme mod- 
els, corresponding to the brightest and the faintest thick 
disks that are compatible with the photometric data. For 
the "faint thick" model, we assume that the thin disk is 
larger and brighter, by decreasing fiojhin an d increasing 
h r ,thin and Zojhin by their uncertainties, and that the 



thick disk is smaller and fainter, by adjusting the same 
parameters in the opposite direction. We then recalcu- 
late fthick- For the "bright thick" model, we make the 
opposite parameter shifts. The "faint thick" and "bright 
thick" models are truly extreme fits to the photomet- 
ric data, and result in reduced chi-squared (x 2 ) values 
that are statistically much worse than the "preferred" 
models (x 2 - preferred > l)i particularly for the "faint 
thick" model. For completeness, we also derive a "sim- 
ple" model assuming that the midplane and above plane 
spectra are entirely dominated by the thin and thick 
disks, respectively. These latter fits are shown in Figure 
H3 We do not include radial gradients in fthick, because 
our photometric decompositions suggest the thick and 
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Fig. 5. — Extracted rotation curves for FGC 1415, showing the rotation along the midplane (left) and above the midplane plane (right). 
The format is the same as Figure 0] There are indications of a slight lag in the off-plane rotation curve. 



thin disks have comparable scale lengths. Likewise, we 
neglect seeing, which should have no effect on the mid- 
plane value of f thick- Seeing affects the off-plane values 
only marginally, shifting the value of f thick by ~ 0.05, 
which is an order of magnitude below the range of f thick 
values we explore. 

We bin our model rotation curves using the same bin 
sizes and flux weighting as our observations. The rota- 
tion curve parameters of both disks are then adjusted to 
minimize x 2 > subject to the constraints that both rota- 
tion curves have the identical central velocity vq, that 7 is 
greater than 0.7, and that the galaxy's kinematic center 
ro is fixed to the flux maximum of the galaxy's contin- 
uum emission. The value of v Ci thick is also constrained 
to be greater than -200 km s _1 to prevent divergent fits. 
We therefore have 7 free parameters: vq, u c ,thin, ft.thim 
7thin, w c ,thick, r t , thick, and 7thick- The best fit values are 
listed in Tables 01 & El along with Xtotali which is the to- 
tal reduced chi-squared for the photometric model plus 
rotation curve fit. The "faint thick" disk models fail to 
converge to reasonable values, but we include the fits for 
completeness. The resulting rotation curve fits are shown 
in Figures El El and, |H1 For both galaxies, every model 
has the rotation of the thick disk lagging the thin disk, 
or even counter-rotating. 

We now test the justification for modeling the rota- 
tion curve data with two distinct kinematic components. 
Specifically, we are attempting to rule out the hypothesis 
that a single underlying rotation curve can explain both 
the midplane and above plane observations. Using \ 2 
as a measure of goodness-of-fit, we expect the minimized 
values of x 2 to be drawn from a chi-squared distribution 
for N — M degrees of freedom, where ./V is the number of 
data points and M is the number of adjustable parame- 
ters in our fit. We calculate the probability Q that a value 
of x 2 larger than the one we find should occur by chance. 
When the midplane and off-midplane are fit with inde- 
pendent rotation curves (the "simple" model), we find 
Q=94% for FGC 1415 and 91% for FGC 227, suggesting 
that this is a very good fit and/or we have overestimated 
our uncertainties. We then compute the probability that 
the off-plane observations are described by same simple 
thin disk rotation curve that fit the midplane, but binned 
in the same manner as the thick disk observations. This 



tests if both the midplane and the off-plane observations 
could be sampling the same underlying rotation curve. 
We find that the off-plane observations of FGC 1415 have 
a 24% chance of matching the thin disk rotation curve. 
For FGC 227, there is only a 4% chance the off-plane ob- 
servations are fit by the simple thin disk rotation curve 4 . 
Although these probabilities do not formally eliminate 
the possibility that the midplane and off-plane obser- 
vations are due to the same underlying rotation curve, 
they do suggest that it is highly unlikely. When coupled 
with the additional facts that our photometric decompo- 
sitions strongly support two disk components, that there 
are different colors of the stellar populations at the two 
slit positions, and that the kinematics of the thick disk 
and thin disk kinematics in the Milky Way are different, 
we believe that our assumption of two kinematic compo- 
nents is strongly justified. 

A final concern is that our observations could be biased 
by the effects of dust in the galaxies. Obscuration by dust 
may systematically reduce the velocities measured near 
the centers of galaxies, where the optical depth due to 
dust is largest. These central velocities will be systemati- 
cally lower due to the failure of observations to penetrate 
to the center of the disk, allowing only stars with small 
line-of-sight velocities to contribute to the spectra. The 
resulting obscured rotation curve would then appear to 
be that of a rotating ring, and would thus mimic solid- 
body rotation a t the inner radii th at are strongly ob- 
scured by dust. iKregel et alJ l)2004l) detect just such a 
signature in longslit spectroscopy of NGC 891, a massive 
galaxy with a prominent dust lane. 

The stronger the extinction caused by dust, the more 
an obser ved rotation curve wou ld approach solid-body 
rotation l)Goad fc Roberts! 11981(1 . This could be cause 
for concern as our RC's are slowly rising and fairly lin- 
ear. Using Monte Carlo radiat ive transfer techn iques, 
iMatthews fc Woodl po"Ol|) and iBaes et all p00^) find 
that realistic optical depths cannot explain the slowly 
rising rotation curves in edge-on LSB galaxies such as 
ours, which lack dustlanes and have a very patchy, highly 

4 These results explain the failure of the "faint thick" disk models 
to converge. The "faint thick" disk model assumes that the thin 
disk dominates the light at both slit positions, which is essentially 
equivalent to the hypothesis we are ruling out. In such a case, the 
thick disk is driven to extreme values in order to minimize x 2 - 
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Fig. 6. — Rotation curve fits to the binned rotation curves for FGC 227 (left) and FGC 1415 (right), assuming that the midplane (solid 
line) and off-plane (dotted line) slit positions sample only the thin and thick disk components, respectively (fthick = in the midplane 
and fthick = 1 above the plane, corresponding to the "simple" model in Tables l3*lfcl4*l. Top panels show the derived rotation curves, while 
bottom panels show those curves binned and flux- weighted identically to our observations. Points in the lower panels are the same as in 
Figures [3] and [S] Dashed lines show the W50/2 velocity for each galaxy. 
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Fig. 7. — Derived rotation curves using our preferred photometric decompositions. For FGC 227, we assume the midplane flux is 20% 
thick disk light while the off-plane observations are 50% thick disk light. For FGC 1415, we assume the midplane flux is 30% thick disk 
light while the off-plane is 70% thick disk light. The lower panels display the rotation curves after they have been binned and flux-weighted 
in the same manner as our observations. Points in the lower panels are the same as in Figures [I] and [F] Dashed lines show the W50 /2 
velocity for each galaxy. 



clumped distribution of dust ijDalcanton et al.1 I2004|h 
IMatthews fc Wood! IpOOlJ) conclude that dust can be ig- 
nored even at incl inations of 90° in LSB galaxies and 
iBosma et al.l l|1992f) find that galaxies with rotational ve- 
locities ~ fOO km s~ x are transparent at optical wave- 
lengths. Our galaxies ar e co mparable to t hose in the 
IMatthews fc Wood! l)2001j) and lBosma et alJ Jl992) stud- 
ies and thus we are confident that our spectra are a true 
measure of stellar orbits and not an artifact of dust ex- 
tinction. 

There are three additional pieces of evidence suggest- 
ing that our results are largely unaffected by dust. First, 
our galaxies are bluer in their midplanes than above 
them, indicating that they do not have central dust lanes. 
Second, our observations show that the off-plane rotation 
curves are lagging the midplanc-the opposite of what we 
would see if there was a dusty midplane. If dust were 
affecting the midplane, then the thick disk would be lag- 
ging the thin disk by even more than we have reported. 
Finally, our observations, while not in the far infrared, 
are red enough (~ 8000 A) that we should be able to 
penetrate the small amount of patchy dust which un- 



doubtedly exists in the galaxies even in the absence of a 
concentrated dustlane. 

5. VELOCITY DISPERSION 

We now turn to an analysis of the line-of-sight velocity 
dispersion (LOSVD) of our galaxy spectra. While di- 
rect fitting returns quality fits for recessional velocities, 
we find that at low S/N this method fails when adding 
additional free components, and is thus unable to mea- 
sure the LOSVD. To measure the LOSVD we have there- 
fore used the updated cross-correlation (XC) method of 
iStatlerl <|1995lh This method measures the LOSVD by fit- 
ting the peak of a template star-galaxy XC with a broad- 
ened stellar autocorrelation function (AXF). 

To avoid possible errors from template mismatch, we 
use two stellar templates, HD4388 (spectral type K3III) 
and HD2I30I4 (spectral type G9III). Both were observed 
with GMOS-N for other programs and were downloaded 
from the Canadian Astronomy Data Centre (CADC) 
archive 5 . The stars were reduced using standard IRAF 
routines and binned logarithmically in wavelength. Spec- 

5 Guest User, Canadian Astronomy Data Centre, which is oper- 
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Fig. 8. — Shaded regions showing the full range of rotation curve fits. The lighter region represents the full range of model thin disks, 
and the dark region shows the range of thick disks, based on Tablcs|3&2] The solid lines represent the rotation curves calculated assuming 
the preferred fthick values based on photometric decomposition. 



TABLE 3 

FGC 227 Rotation Curve Parameters 



Model 


midplane, fthick 


off-plane,/ tWcfc 


v 


^"t,thin 


^t,thick 


Ithin 


"fthick 


^c,thin 


"e, thick 


^total 


preferred 


0.2 


0.5 


5493. 


4.6 


0.9 


20. 


1.4 


in. 


-51. 


5.0 


faint thick 


0.1 


0.1 


5493. 


4.8 


0.5 


200. 


207.5 


87. 


-37. 


15.8 


bright thick 


0.5 


0.9 


5493. 


5.0 


5.8 


28. 


10.6 


145. 


16. 


8.2 


simple 


0.0 


1.0 


5493. 


5.1 


7.0 


176. 


12.0 


80. 


35. 





tra of reference arcs were used to determine the amount 
of artificial broadening needed to match the resolution of 
the stellar template to the galaxy spectrum. 

Before performing the XC, we performed several pre- 
processing steps. First, we once again extracted galaxy 
spectra by summing our reduced 2D image along the 
spatial direction until adequate S/N was reached (~ 10). 
Next, the galaxy and stellar templates were normalized 
by dividing by a low-order fit to the continuum, af- 
ter which the mean was subtracted. The spectra were 
padded (100 pixels) to eliminate "wrap-around" effects 
when cross-correlating. Residual low frequency compo- 
nents were removed by tapering the ends of the spectra 
in the Fourier domain. 

The faintness of our off-plane spectra requires that we 
bin or observations spatially to reach adequate signal- 
to-noise. However, using such large extraction bins ar- 
tificially broadens our galaxy spectra by including stars 
from different parts of the rotation curve. To eliminate 
this, we first run the XC fitting code to find the reces- 
sional velocity along the galaxy. We then shift individual 
rows of the 2D galaxy spectrum to remove the measured 
rotation and re-extract ID spectra. The XC fit is then 
run again in order to measure the LOSVD. We find no 
major effects from template mismatch and can use our 
template stars interchangeably. When fitting for rota- 
tion, the two templates return nearly identical results 
with an RMS spread of 1.2 km s" 1 with a maximum of 
4 km s _1 difference. Velocity dispersions show a ~ 5 km 
s _1 systematic difference which is less than the statisti- 
cal uncertainties and may result from mismatching the 

ated by the Herzberg Institute of Astrophysics, National Research 
Council of Canada. 



dispersions based on reference arc spectra. 

We must emphasize that the velocity dispersion mea- 
surements are very uncertain, especially for FGC 227. In 
the direct fitting method, we were able to mask trouble- 
some skylines and use partial sections of the Ca triplet. 
When using the XC method, masking small sections of 
spectra can create large changes in the Fourier domain, 
so we were forced to use only 2 of the 3 Ca absorption 
features for FGC 227. Therefore, the LOSVD measure- 
ments for FGC 227 should be viewed with caution. 

The resulting measurements of the line-of-sight veloc- 
ity and velocity dispersion are shown in Figure El The 
rotation curves from this method confirm the results of 
our direct fitting method. Both methods show a possi- 
ble slight lag for the thick disk of FGC 1415 and zero 
net rotation for the thick disk of FGC 227. The LOSVD 
measurements of FGC 1415 are similar for the midplane 
and off-plane, and show a roughly constant velocity dis- 
persion with radius. The midplane of FGC 227 shows a 
possible increase in the LOSVD with radius. Unfortu- 
nately, the off-plane observations of FGC 227 have suffi- 
ciently low S /N and high sky-line contamination that the 
LOSVD measurements have very large error bars. We are 
therefore unable to detect the increase in velocity disper- 
sion expected for a thick disk which is not supported by 
rotation. The midplane LOSVD is qualitatively consis- 
tent with the presence of a counter rotating thick disk, 
but the off-plane observations have too low S /N to make 
firm conclusions. 

To understand the radial behavior of the LOSVD, we 
modeled the velocity fields of two superimposed rotat- 
ing disks viewed edge-on, corresponding to a thick and 
thin disk rotating at two different speeds. We have as- 
sumed each disk has an intrinsic velocity dispersion of 
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TABLE 4 

FGC 1415 Rotation Curve Parameters 



Model 


midplane, fthick 


off-plane, fthick 


vo 


ft, thin 


ft,thick 


'Jthin 


'thick 


^c,thin 


Vc, thick 


^total 


preferred 


0.3 


0.7 


1522. 


2.8 


0.8 


17. 


0.7 


102. 


28. 


24 


faint thick 


0.2 


0.4 


1522. 


2.7 


21.4 


100. 


15.2 


97. 


-200. 


45 


bright thick 


0.4 


0.9 


1522. 


2.8 


2.5 


21. 


5.8 


101. 


40. 


26 


simple 


0.0 


1.0 


1522. 


2.7 


2.2 


33. 


13.4 


77. 


41. 
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Fig. 9. — Results from measuring line of sight velocity and velocity dispersions using the improved cross-correlation method. Dashed 
lines in the LOSVD plots show our pixel widths. Dashed lines in the velocity plots are shown for reference between midplane and off- 
plane components. Any dispersion below that level would be poorly sampled and we would expect the XC method to fail. The rotation 
measurements are consistent with our least-squares method while FGC 227 shows a suggestive increase in LOSVD with radius. 



~ 30 km s _1 . We created a series of artificial spectra 
from a combination of two broadened template stellar 
spectra with different flux strengths and recessional ve- 
locities, and then recovered their LOSVD with our XC 
procedure. A sample of the resulting LOSVD are plot- 
ted in Figure ^| We find that, at our resolution, the 
XC method is unable to separate the cross-correlation 
peaks unless the two stellar components have velocity 
differences > 200 km s _1 . For a case like FGC 1415, 
where the thick disk is only mildly lagging the thin disk, 
we would therefore expect only a slight increase in the 
observed LOSVD to occur with radius, consistent with 
the constant LOSVD seen in Figure [5] However, we find 
that adding a second component at a different velocity 
can dramatically increase the observed LOSVD if the ro- 
tational velocity difference is in the range 90-150 km s _1 
and the second component comprises more than 20% of 
the total flux. For FGC 227 we observe a large increase 
in velocity dispersion with radius as expected if there is 
a counter- rotating thick disk present. 

6. DISCUSSION 

The results of 21 and [5] indicate that the thick disks of 
FGC 227 and FGC 1415 have very different kinematics. 
The raw rotation curves in Figures show that both 

galaxies appear to be rotating more slowly above the 
midplane than within it. However, when the thick and 
thin disk kinematics are decomposed (Figure |HJ, even 
stronger differences emerge. 

In FGC 1415, the decomposition strongly favors a thick 
disk that rotates in the same direction as the embedded 
thin disk, but with a lag of 50% or more. The kinematics 
of FGC 1415 are therefore fully consistent with those 
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Fig. 10. — The expected observational effect of two disk compo- 
nents with different rotational velocities on the observed LOSVD. 
Plot was made by combining two artificial spectra with intrinsic 
velocity dispersions of 30 km s — 1 . For galaxies like FGC 1415 
with mildly lagging thick components (max lag ~ 80 km s — 1 ), 
we would expect a LOSVD enhancement of ~ 15 km s — 1 in both 
the midplane and off-plane observations. If FGC 227 has a counter 
rotating component we would expect a larger increase in LOSVD 
with radius. 



seen in the Milky Way, where the thick disk lags the 
rotation of the thin disk by up to 100 kms" 1 at lar ge scale 
heights (|Oilmore et al J 12003 lParker"eTaIll200l . This 
suggests that any of the theories that can explain the 
kinematics of the Milky Way's thick disk can also explain 
the kinematics of FGC 1415. 

The kinematics of FGC 227 are quite different. The 
very slow rotation above the plane in FGC 227 is due 
entirely to contamination from the thin disk. Unlike 
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FGC 1415, the thick disk in FGC 227 is counter-rotating, 
if it is even rotating at all. While monolithic collapse sce- 
narios could potentially explain the rotation of FGC 1415 
and the Milky Way's thick disk, they are unable to ex- 
plain the slow counter-rotation seen in FGC 227. In such 
scenarios, the rotational lag between the thick and thin 
disks arises from the spin up of gas as it collapses into 
the thin disk while preserving its angular momentum. 
Thus, the thick and thin disks must always rotate in the 
same direction, in contrast to the behavior of FGC 227. 
The data for FGC 227 could marginally accommodate a 
very slowly co-rotating thick disk, but only if its thick 
disk formed very early in the collapse process, before the 
gas had spun up substantially. This would form a thick 
disk that is substantially thicker and more extended than 
the thin disk. However, the structure of FGC 227's and 
FGC 1415's thick and thin disks are comparable, ruling 
out monolithic collapse as a possible explanation. 

The data for FGC 227 are also incompatible with 
models where the thick disk forms via vertical heat- 
ing of a young thi n disk. Simulations show that 
infalli ng satellites llQuinn et al.l 1 19931 I Walker et al.l 
Tflflfit iVelazquez fc White! 119991 lAguerri et al] 120011: 
Benson et al.ll2004D and dark mat ter halo substructure 
((Font et al.ll2001t lArdi et al.ll2003j) increase the velocity 
dispersion of the disk stars, but do not substantially al- 
ter their angular momenta. The slowly rotating thick 
disk of FGC 227 is therefore unlikely to have formed 
from a rapidly rotating thin disk. Another possibil- 
ity is that the thick disk formed from a thin disk that 
was also slowly rotating. Such a disk must have had 
a much smaller scale length because of its low angu- 
lar momentum, and thus, the heated thick disk would 
have a smaller scale length as well. This is in con- 
trast to observatio ns that consistently find large thick 
disk scale lengths flWu et all l2002t INeeser et aj J 12001 
lYoachim fc Dalcantonl 120051 iPohlen et alJ 12004(1 We 
therefore rule out vertical heating as a plausible mecha- 
nism for forming FGC 227's t hick disk. Ou r concl usions 
are consistent with those of IPohlen et alJ 1(20041) . who 
performed a very different analysis of thick disk struc- 
tural parameters in SO galaxies. 

Having ruled out monolithic collapse and vertical heat- 
ing, the only remaining possibility is that the thick 
disk formed from direct accretion of stars from satel- 
lites. In this scenario, two or more galactic fragments 
merge together to form the final galaxy. The majority 
of the stars in the satellites remain in the thick disk, 
due to their inability to lose kinetic energy and "cool" 
into a thinner disk. Meanwhile, the gas in the merg- 
ing satellites collapses further into a rotating disk, con- 
verts into stars, and forms the thin disk. The gas in the 
thin disk may be supplemented at later times by addi- 
tional gas infall. These "shredded satellites" are becom- 
ing popular as e xplanations for a wide variety of struc- 
tures in the MW llNavarro et alJ200llMartin et al.l2004t 
iParker et aDl2004HHelmi et allbomli . as well as for low 
veloci ty stars at hig h latitude in the Mil ky Wa y thick 
disk ((Gilmore et al.ll200l iSoubiran elaD 12003]) . This 
class of models can also reproduce the enrichment pat- 
terns seen in Milky Way halo stars ((CI ribA^B£exj|2X)0i]|; 
Nisser, fc Schuste1ll997t iBTook et a,Ul200.3t ISa.leh et alJ 
2004). 

Recently, specific satellite accretion models for the for- 



mation of the thick disk have been advocated on the 
basis of numerical simulations. lAbadi et alJ l|2003D and 
iBrook et all (|2004D both identify thick disk components 
in gala xies formed in cosm ological simulations. The thick 
disk in lAbadi et alJ l(2003l) is composed of s tars accreted 
from merging galactic "sub-units" , while IBrook et alJ 
(2004) find that the thick disk stars form during the 
chaotic merging of gas-rich clumps that produce the disk 
itself. These new simulations form thick disks through 
a series of mergers, in contrast to older models that re- 
lied on a single fi nal merging event to heat a thin disk 
into a thick disk iQuinn et al.lll993T: iWalker et"al]ll99ft 
IVelaznuez White! l1999f Th e merging scenario ex- 
plored bv the lAbadi et all (|2003) and lBrook et all l|2004ft 
simulations may sometimes reproduce the behavior of the 
older disk heating models, for example, when one mas- 
sive satellite comes in at late times. However, the newer 
models provide mechanisms for producing a wider range 
of thick disk properties. Variations in the orbital proper- 
ties and gas richness of the merging proto-galactic clumps 
easily yield variations in the kinematics of the thick disk, 
as well as its size and its luminosity relative to the thin 
disk. 

For example, satellite merging can readily produce a 
modest rotational lag for the thick disk. In general, the 
same satellites produce both the gas that settles into the 
thin disk and the stars which are left in the thick disk, 
assuming little subsequent gas accretion. The thin disk 
and thick disk should therefore show similar angular mo- 
mentum distributions, but the rotation speed of the gas 
should be slightly higher, because it has contracted fur- 
ther into the halo. This process naturally leads to a thin 
disk which rotates more rapidly and has a smaller scale 
length than the thick disk, as seen in the Milky Way and 
FGC 1415. 

A different combination of orbital parameters and gas 
richnesses in the merging satellites can also produce the 
slowly rotating thick disk of FGC 227. If the merging 
satellites had little net angular momentum (for exam- 
ple, two equal mass satellites orbiting the center of mass 
in opposite directions), then the stars the satellites de- 
posited in the thick disk would show little net rotation. 
However, if one of the merging satellites was particularly 
gas rich, then it alone could provide enough gas to cre- 
ate the rapidly rotating thin disk. This pathway would 
lead to a rapidly rotating thin disk, a stationary or even 
counter-rotating thick disk, and comparable scale lengths 
for the thick and thin disks, exactly as seen in FGC 227. 
The resulting thick disk would also have a large velocity 
dispersion, but unfortunately, our data are not sufficient 
to conclusively demonstrate this for FGC 227. 

One possible limitation with the stochastic satellite ac- 
cretion model is its questionable ability to consistently 
form highly flattened stellar systems. Our past imag- 
ing indicates that late-type disks rout inely host flat- 
tened thick disks with ~5:1 axial rat ios l | Morrison et all 
[ 19971 iDalcanton & Bernstein! l200l INeeser et alJ 120021: 
IWu et afl l2002'l. However, if satellites are accreted from 
random directions, then there is no a priori reason why 
the the satellite debris should necessarily be aligned with 
the embedded thin disk. One would therefore expect to 
see a significant fraction of galaxies that lack thick disks, 
but are instead embedded in round stellar halos. 

There are several possible resolutions to this problem. 
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The first is that much of the satellite debris that is not 
aligned with the disk instead becomes identified with the 
bulge. This may be particularly true when the merging 
satellites are gas-rich and are on largely radial orbits, 
leading to centralized bursts of star formation and the 
formation of a high surface brightness bulge. Any stars 
formed during the burst, or accreted later, will tend to be 
identified with the brighter bulge component. However, 
this potential solution cannot explain why flattened thick 
disks are seen in bulgeless late-type galaxies. 

The second possibility is that the orientation of the an- 
gular momentum vector of the accreted thick disk stars 
may be preferentially aligned with that of the accreted 
gas that forms the thin disk, even if the amplitude of the 
angular momentum vectors differ. lAbadi et al.l (2003) 
have argued that this alignment naturally arises from the 
tendency of satellites to travel along filaments, providing 
a preferred direction for satellite accretion. The align- 
ment also may occur through dynamical friction. Satel- 
lites on circular orbits will have longer lifetimes, and will 
thus be in fluenced more by d ynamical friction. Simu- 
lations by W alker et alJ l)199 6i show how a non-planar 
satellite orbit can be quickly dragged into the plane by 
dynamical friction in only ~ 1.5 orbits. The relative con- 
tri bution of stars from satellites with different orbits in 
the lAbadi et alJ (|2Q03) simulations support this picture. 

The third possibility is that the tendency to detect 
flattened thick disks is a selection bias against detect- 
ing rounder distributions of satellite debris. If a merg- 
ing satellite does not happen to be in a coplanar orbit 
with the thin disk, then it will deposit stars over a larger 
volume throughout the stellar halo. The resulting dis- 
tribution of stars will have a much lower surface bright- 
ness than if the accreted stars were concentrated near 
the midplane in a thick disk. This difference naturally 
biases against detecting a rounder distribution of stellar 
debris. IWu et all l)2002|) may have detected this compo- 
nent in deep intermediate-band imaging of NGC 4565, 
as may have iZibetti et all l)2004f) using stacked images 
of edge-on galaxies from the Sloan Digital Sky Survey. 
However, there is a possibility that some of this light is 
due to the extended light from bulges embedded in the 
coadded thin disks, so this latter result remains ambigu- 
ous. A final possibility is that the surface brightness of 
any round component will be further reduced by orbital 
instabilities. Assuming that the potential of most galax- 
ies are triaxial, orbits around the intermediate axes will 
be unstable. These orbits will tend to reorient around 
the stable short or long axes, pulling a fraction of the 
stars back into a coplanar orientation. 

7. CONCLUSION 

We have measured the rotation curves and velocity dis- 
persion of stars at and above the midplane of two edge-on 
galaxies. Comparisons between the two positions indi- 
cate that the stars above the plane are significantly lag- 
ging the rotation of the stars in the midplane. We have 
decomposed the observed rotation into two components, 
corresponding to the thick and thin disks. We find that 
the rotation speed of the thick disk component is less 
than 50% of the rotation speed of the thin disk, and is 
counter-rotating in one of the two cases. We use these 
observations to conclusively rule out older monolithic col- 
lapse or disk heating models for forming all thick disks. 



We conclude that the only viable mechanism capable of 
forming counter-rotating thick disks is one where the ma- 
jority of thick disk stars are directly accreted from merg- 
ing satellites. Simulations suggest that the early merging 
required by these models is a generic feature of hierar- 
chical galaxy formation. As such, the models imply that 
all disk galaxies should host thick disks. This suggests 
that merger models for the thick disk not only account 
for the variation in kinematics, but also explains the ap- 
parent ubiqui ty of thick disks around disk galaxies, as 
argued for bv lDalcanton k, B crnstciiJ l)2002j) . 

We are currently expanding the sample to encompass a 
wider range of galaxy masses. We will discuss the larger 
sample and the measured velocity dispersions in a future 
paper. 
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